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Abstract 

Using the NA48 detector at the CERN SPS, 31 Ks — > 7r°77 candidates with an estimated background 
of 13.7±3.2 events have been observed. This first observation leads to a branching ratio of BR(Ks — > 
7r°77) = (4.9 ± l.Gstat ± 0.9 sy st) x 10~ 8 in agreement with Chiral Perturbation theory predictions. 
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1 Introduction 

Radiative non-leptonic rare kaon decays have proved to be useful tests for low energy hadron 
dynamics studied in the framework of the Standard Model by the Chiral Perturbation theory (xPT). 
This paper describes the first observation of the K$ — ► tt°jj decay obtained from data taken by the 
NA48 experiment in the year 2000. 

In the decay Ks — » 7r°77, the photon pair production is described by an amplitude dominated by 
a pseudo-scalar meson pole. In xPT this pole is dominated by ir° contribution, and the lowest order 
amplitude is non- vanishing, in contrast to the similar Kl — * 77 decay. The theoretical prediction for the 
branching ratio is 3.8 x 10~ 8 with higher order corrections expected to be small [T] and is quoted in the 
kinematic region z — n&„/m K > 0.2 which is free from the overwhelming Ks — > ir°ir° background. 
A measurement of the branching ratio can provide information about the presence of extra non-pole 
contributions studied e.g. in 2 . In addition, the momentum dependence of the weak vertex which is 
predicted by xPT can be tested by the measured shape of the z spectrum. The lowest previously published 
limit on the branching ratio is BR(Ks — ► 7]"°77)z>o.2 < 3.3 x 10~ 7 at 90% confidence level [3]. 

2 Experimental set-up and data taking 

The NA48 detector was designed to measure direct CP violation in the decays of Kl and Ks into 
7T7T, described by the parameter Re(e'/s), using simultaneous far- and near-target beams The analysis 
presented in this paper is based on data recorded during a special 40-day run performed in the year 2000 
with near-target beam only and with an intensity of about ~ 10 10 400 GeV protons hitting the target 
during the 3.2 s long SPS spill. The exit of the collimator, 6 m downstream of the target, is followed 
by a wide vacuum tube approximately 113 m long containing the decay region and terminated by an 
aluminium window close to the detector 1 ' . 

The detector elements used in the analysis are the following: 

— A liquid krypton calorimeter (LKr) 0, placed less than 2 m behind the aluminium window, is used 
to measure the energy, position and time of electro-magnetic showers. The energy resolution is 

a{E)/E ~ O.O9O/£0O.O32/\/£; 0.0042, (1) 

where E is in GeV. The position and time resolutions for a single photon with energy larger than 
20 GeV are better than 1.3 mm and 300 ps, respectively. 

— A sampling hadron calorimeter composed of 48 steel plates, interleaved with scintillator planes, 
with a readout in horizontal and vertical projections. 

— Seven ring shaped scintillator counters equipped with iron converters (AKL), used to detect photons 
escaping the outer limits of the calorimeter acceptance. 

The trigger decision, common to 7r°77 and 7r°7r n decays, is based on quantities which are derived from 
the projections of the energy deposited in the electro-magnetic liquid krypton calorimeter UJ. The trigger 
required that the total deposited energy E to t is larger than 50 GeV, the radial distance of the centre of 
energy from the beam axis is smaller than 15 cm and the proper life time of the kaon is less than 9 Ks 
lifetimes downstream of the collimator. The inefficiency of the trigger was measured to be at the level of 
0.1%. 

More about the detector and the experimental configuration during data taking in the year 2000 
can be found in [7|. 

3 Event selection 

The energies and positions of electro-magnetic showers initiated by photons, measured in the liquid 
krypton calorimeter, are used to calculate the kaon energy and decay vertex. To select Ks — > ir°jj and 
Ks — > 7r°7r° candidates, all events with > 4 energy clusters are considered. All combinations of four 
clusters which are within 5 ns from the average time and with no other cluster with energy > 1.5 GeV 
closer in time than 7 ns with respect to the event time are selected. The event time is computed from the 
times of the most energetic cells of the selected clusters. In addition, the clusters must pass the following 
cuts: 

— The energy of each cluster must be greater than 3 GeV and less than 100 GeV. 

— The transverse distance between two clusters is required to be greater than 10 cm. 

— The total energy of the selected cluster combination is required to be less than 170 GeV and greater 
than 70 GeV. 

x ) In the usual NA48 experimental configuration the vacuum tube is terminated at 89 m by a thin Kevlar window followed 
by a helium filled tank which contains drift chambers and a spectrometer magnet. For these data however, the Kevlar 
window and the drift chambers were removed, the helium tank was evacuated and the magnet was switched off. 
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— The distance of the centre of energy from the beam axis, 



p V(E.^) 2 + (E t ^) 2 , 9 v 
Rc = E^ (2) 

is required to be less than 5 cm, where Ei, Xi, yi are the i-th cluster energy, x and y coordinates at 
LKr, respectively. 

— The energy deposited in the hadron calorimeter must not exceed 3 GeV in a time window of ±15 
ns around the event time. 

— The AKL counter should not register any hit in a time window of ±7 ns around the event time 
The decay vertex position z ver t e x of a kaon is calculated using the kaon mass constraint, rnx, 



■^vertex — ^Lkr \ t ->) 

m K 

where ZLkr is the longitudinal coordinate of the LKr calorimeter with respect to the end of the 
collimator and dij is the transverse distance between i-th and j-th cluster at the calorimeter. The invariant 
mass of photon pairs is calculated using z ver tex- 

The Ks — * 7r°77 candidates must have a photon pair with an invariant mass within 2 MeV/c 2 of 
the 7T° mass, m^o, and the other pair must satisfy the condition mL^/mf^- > 0.2. In the K$ — > 7r°7r° 
sample a \ 2 cu t of 27 (~ 5<j) is applied to the invariant masses of photon pairs, defined as 



x 2 



(mi + m 2 )/2 - m K o 



(mi - m 2 )/2 



(4) 



where a± are the resolutions on (mi ± m-i)l1 measured from the data and parametrised as a function of 
the lowest photon energy. The typical value of a + is 0.4 MeV/c 2 and cr_ is 0.8 MeV/c 2 0. For Ks — * ir°jj 
candidates the \ 2 is required to be larger than 5400. 

In order to suppress the large background from Ks — * 7r°7r2> decays with one particle escaping 
the acceptance the decay region was restricted to be between -1 m and 8 m with respect to the exit 
of the collimator. Most of the Ks — * ^ a ^% background events have an apparent z vertex downstream of 
this region because of the missing energy. This condition is even more effective against background from 
Kl — > 7r 7r°7T . The Ks — > 7r° 7r^ background is further suppressed by imposing the 7r° mass hypothesis 
on the remaining four pairings 2 ) of photon showers and calculating the vertex positions z^o for each of 
these 7 pairs. Denoting as z* the closest z^o to the exit of the collimator, events with 

z* < -4.5m or z* > z vertex + 10m (5) 

are accepted (Fig.^l. This cut is very effective in suppressing Ks — > ^ 0/ ^% background because, as shown 
by Monte Carlo studies, this background is completely dominated by events with one of the it° D — > eej 
decay products escaping detection. In addition, this cut removes events from inelastic scattering of beam 
particles in the collimator with multiple ir° production. The downstream end of this cut is intentionally 
extended in order to reduce background from pile-up of a ir° from Ks — > 7r°7r° with 7's from another 
decay. 

A small amount of Ks — * it tTj^ background survives the above cuts if a 7 from one ir° and a 7 or 
an e from the other tt° overlap. In this case the reconstructed z ver tex does not move downstream because 
the entire kaon energy is collected in the LKr calorimeter and there is no pair of showers which would give 
a correct z n o under a ir° hypothesis. Still, by assuming this type of overlap the correct vertex position 
can be reconstructed by 



Zoverlap ^-LATr 



,2 E, I ./ ( E i E J d ij 



fll<*14 3 ( 6 ) 



-\ 14 Ef^i 



where index 1 refers to the 7 from the 7r°, indexes 2 and 3 to 7 or e's from tt d and index 4 to the 
overlapping shower. For 7r°77 events z over i ap ^ z vertex and hence a cut \z over i ap — z vertex \ > 1.5 m on all 
12 possible combinations reduces this type of background to a sufficiently low level, without significant 
loss of signal. 

2 ' If indexes 1 and 2 refer to 7's assigned to the n° then these four pairings are 1-3, 1-4, 2-3 and 2-4. 
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Figure 1: Illustration of z vertex and z* cuts on the data. All other cuts were applied. Rejected areas are 
shaded. Both variables are defined with respect to the exit of the collimator. The accumulation of events 
around small z* values is due to K$ — > 7r ° 7r i) background at high z vertex and due mainly to pile-up 
events at low or negative z ver t e x- 



In order to remove events with hadrons or overlapping showers from the 7r°77 sample, the shower 
width is required to be less than 3er above the average value for photon showers of a given energy. This 
cut, which is calibrated using showers in Ks — > 7r°7r° decays, removes < 1% of good K$ — > ir°jj events. 

In addition, residual background from S° — * Air with subsequent A — * nir° decay, where one 7 
escapes and the neutron produces a photon- like shower in the LKr calorimeter, is suppressed by exploiting 
the large momentum asymmetries in both S° and A decays and accepting only events with: 

\E 3 - EA , (E 3 + EA - (Ei + E 2 ) 

1 J 1 < 0.35 and V ^ if 1 ' < 0.3 (7) 



E 3 + £4 Ei + Ei + Ej, + E4 

where indexes 1 and 2 refer to the shower pair identified as photons from the 7r° and 3, 4 to other two 
showers. 

The background from accidental pile-up events is reduced by strengthening the requirements on 
shower times. Two configurations of a pile-up are considered: 3+1 showers and 2+2 showers. The 3+1 
configuration is best described by a variable 

max (8) 

where the difference between any shower time and the average of the other three is maximised. Similarly 
for the 2+2 configuration one can use a variable 

^2m,ax — [pi ^ ^ tj/^]max (9) 

in which only two of three remaining shower times are averaged. For the 3+1 pile- up topology ti max 
and t 3max are always equal within the time resolution, however in the 2+2 case the two variables have 
different values in principle. This allows one to analyse this type of background in more detail. For the 
selection of the signal sample a single cut on i 2mal of <1 ns is chosen, because t^max describes both 
configurations in the same way. 

After all selection cuts 31 events remain in the sample. For the normalisation a sample scaled-down 
by factor 1000 of about 285 x 10 3 K$ — > 7r°7r events has been used. 



4 Background subtraction 

The following processes have been considered as potential sources of residual background: 

— irreducible Kl ~* ir°j~f background from the Kl component of the beam 

— Ks —>■ 7r°7r° with mis-reconstructed energy 

— Ks — > 7r°7T£) with a 7 or an e escaping detection or with a shower overlap 
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— Hadron background, mainly S° — > Air with subsequent A — * nir° , with three photon showers and 
one narrow neutron shower in the LKr calorimeter 

— Pile-up of two decays where the two decaying particles originate at the target either from two 
protons (accidental pile-up) or from a single proton (in-time pile-up). 

The amount of Kl - > ir°jj admixture in the signal sample is estimated using the Kl flux measured 
from the Ks — > 7r°7r° rate assuming equal production of Ks and Kl at the target. The acceptance was 
calculated using a Monte Carlo simulation. This background amounts to 3.8 ± 0.2 events. 

Ks — > ir°ir° events can pass the cuts on invariant masses, and especially \ 2 > 5400, only if far 
non-Gaussian tails are present in the energy reconstruction from the LKr calorimeter. In order to study 
cases where one of the four shower energies may be mis-reconstructed, making use of the over-constraint 
of K s -> ir°ir° events, the invariant mass of the event is reconstructed by using only three out of four 
shower energies and all four positions. None of the signal events has been found to have an invariant mass 
close to the kaon mass in any of the shower combinations. In addition, a toy Monte Carlo was employed 
to generate Ks — > 7r°7r° events using realistic resolutions and non-Gaussian tails in energy, enhanced 
by an order of magnitude with respect to those known from E/p studies of electrons from the Re(e' /e) 
analysis. In a sample equivalent to twice the flux of the 2000 near-target run no event passes the x 2 , m -yj 
and z V ertex cuts at the same time. 

The residual Ks — > background has been studied using Monte Carlo simulations. Generating 

1.7 times the collected flux, 4 events have been found to pass all cuts. Three of them have overlapping 
showers and reside in the tail of the z over i ap — z vertex distribution. The fourth event decays in the colli- 
mator, and one of the photons undergoes a conversion in the collimator walls, while the e + e~ pair from 
the 7T° passes trough the central hole of the detector. Scaling these events to the observed kaon flux, 
a background estimate of 2.4 ± 1.2 events is obtained with an uncertainty dominated by the statistics. 
The systematic error of this estimate has been checked by relaxing z* and z over i ap cuts. The number of 
events and their z vertex distribution agree with data within few percent. 

The amount of hadron background surviving cuts on energy deposited in the hadron calorimeter, the 
shower width and the energy asymmetry (J7J was estimated by extrapolating the shower width distribution 
from large widths to the signal region. The shower width distribution of neutrons was extracted from 
fully reconstructed S° — > Air — > nir a ir° events. An estimate of 0.1 ± 0.1 events has been obtained. 

The accidental pile-up background has been studied using time variables timax @ and t^ max (JSJ 
in a time window of 6 ns which is six times larger than the signal time window. In this time window the 
event distribution is not affected by any selection cut. It turns out that about 50% of the background is 
in a 3+1 configuration while 50% is 2+2. Since both of these configurations are described equivalently by 
t2max this variable is used to extrapolate from the control region 2 < \t2max\ < 6 ns to the signal region 
|i2ma:r| < 1 ns assuming a flat distribution. This extrapolation gives an estimate of 7.0 ± 1.3 events for 
the accidental pile-up background. The order of magnitude of this estimate was confirmed by overlaying 
3-photon events from a Ks — > 7r°7r° toy Monte Carlo with random showers detected close in time (20-30 
ns) to fully reconstructed Ks — > 7r°7r events and by taking into account the kaon flux. 



Figure 2: The Rc distribution of data compared to the sum of background models and signal Monte 
Carlo scaled by the calculated branching ratio. The control region for pile-up background subtraction is 
indicated. 




R c (cm) 
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The in-time pile-up background in general has a topology similar to the accidental pile-up, but, 
since it is generated by the same proton in the target, its rate cannot be measured by extrapolating from 
the out-of-time sample. The relative rate of in-time and accidental pile-up has been studied using fully 
reconstructed Ks — * 7r°7r° with additional showers in the LKr calorimeter. It has been found that a 20% 
enhancement of pile-up background due to an in-time component cannot be excluded. In order to estimate 
the amount of total pile-up background without using time variables, the distribution of Rc has been 
employed. The background is extrapolated from a control region 7 < Rc < 10 cm to the signal region 
Rc < 5 cm by using the Rc shape obtained from the out-of-time sample in the t2max control region. This 
extrapolation leads to an estimate of 6.8 ±2.9 events for both in-time and accidental pile-up backgrounds. 
The uncertainty takes into account the statistical error of the data in the control region as well as the 
uncertainty of the extrapolation. This measurement is in good agreement with previous considerations 
and, since it contains the least number of assumptions, it is used for the branching ratio calculation. As 
can be seen from Fig. [2 the Ks — * ^ir^ background is not double-counted in this subtraction because 
it populates the Rc < 7 cm region. 

Recently, the Ks — > 7r°e + e~ decay has been observed for the first time. Using the measured branch- 
ing ratio [H] and taking into account the cutoff on the invariant mass of the e + e~ pair rn\ e /m 2 K > 0.2 
this background has been estimated to contribute 0.6 ± 0.3 events. 

5 Result 

Of the selected 31 events, 13.7 ± 3.2 are estimated to be background (Table [IJ. The probability to 
observe 31 or more events in absence of a signal with a background of 13.7 ± 3.2 events is 1.5 x 10~ 3 . 
Subtraction of the background from the selected sample leads to a signal of 17.3 ± 6.4 events. In the 
control region, defined as 3 < |mi2 — m^ol < 5 MeV, the observed data agree well with the background 
estimate (Fig. |3J). 



Events in 


signal 
region 


control 
region 


Kl — > 7r u 77 background 
Ks — > 7r°7i£j background 
Hadronic background 
Pile-up background 
Ks —> ir°e + e~ background 


3.8 ± 0.2 
2.4 ± 1.2 
0.1 ± 0.1 
6.8 ± 2.9 
0.6 ± 0.3 


0.1 ± 0.0 

4.2 ± 1.6 
0.1 ± 0.1 

4.3 ± 2.0 
« 0.1 


Total background 


13.7 ± 3.2 


8.7 ± 2.6 


Data 


31 
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Table 1: Summary of background composition in both signal and control regions, the latter being defined 
as 3 < |toi2 — m^o | < 5 MeV, compared to data. 




0.132 0.134 0.136 0.138 



0.14 

m 12 (GeV/c 2 ) 



Figure 3: The invariant mass distribution of the photon pair assigned to the tt° . The data in the control 
region are compatible with background while the signal region contains a conspicuous excess indicating 
a Ks — > 7r°77 signal. 
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The detector acceptances of the K$ — > ir°jj decay and of the normalisation decay channel Kg —> 
tt tt q were calculated using a full Monte Carlo simulation of the detector based on GEANT ^U] 3 ). The 
Ks — ► 7r°7r° acceptance was determined to be 18.6%. The K$ — > tt°jj generator used the decay matrix 
element calculated by F using %PT. The acceptance of the Ks — * ir°jj for z > 0.2 is 7.2 ± 0.3%, 
where the uncertainty is extracted from the comparison with a simulation using a pure phase-space decay 
generator. The Ks — > ir°jj acceptance is smaller than that of Ks — ► tt°tt mainly due to cuts on z* Q and 
energy asymmetry defined in (JSJ and Q. 

Normalising the signal to the collected Ks — J ; J ; 
results in a ratio of decay widths 



7r 7r sample and taking into account the acceptances 



T(K S -> 7r°77) 2> o. 2 
T{K S -> 7r°7r°) 



= (1.57±0.51 stat ±0.29 syst ) x 10 
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(10) 



where the systematic uncertainty takes into account the background subtraction and acceptance calcula- 
tion uncertainties. This result is stable against variations of most relevant cuts. In particular, the result 
does not change significantly by releasing the cut on the AKL anti-counter which increases the pile-up 
background by a factor 5, and when requiring no associated hit in a scintillator hodoscope placed upstream 
of the LKr calorimeter, which reduces the Ks — > tt ^ background by a factor 3. Using BR(Ks — > 7r°7r°) 
from [TT], one obtains 



BR{K S 



77°77)2>0.2 



= (4.9 ± 1.6 atot ± 0.9 S j, s t) x 10" 
= (4.9 ± 1.8) x 10" 8 



(11) 
(12) 



which agrees with the predictions of pQ. 

In order to test the momentum dependence of the weak vertex predicted by the xPT the z dis- 
tribution of the sample after background subtraction has been compared to z distributions of simulated 
Ks — > 7r°77 events using two decay generators, one with the xPT matrix element and one with pure 
phase space. Fig. 0] shows that both calculations agree within errors with the data and more statistics 
would be needed to prove the chiral structure of the weak vertex. 
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Figure 4: The z distribution of data after background subtraction compared to a Monte Carlo Ks — > 7r°77 
simulation using phase space (dashed line) and X-P^ matrix element pQ (continuous line) in the decay 
generator. 
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3 ' In order to speed up the simulation, electromagnetic showers are pre-generated and stored in a shower library. 
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